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Abstract 

Wc study Casimir scaling and renormalization properties of Polyakov loops in different ir- 
reducible representations in SU(Af) gauge theories; in particular, we investigate the approach 
to the large- -/V limit, by performing lattice simulations of Yang- Mills theories with an increas- 
ing number of colors, from 2 to 6. We consider the twelve lowest irreducible representations for 
each gauge group, and find strong numerical evidence for nearly perfect Casimir scaling of the 
bare Polyakov loops in the deconfined phase. Then we discuss the temperature dependence 
of renormalized loops, which is found to be qualitatively and quantitatively very similar for 
the various gauge groups. In particular, close to the deconfinement transition, the renormal- 
ized Polyakov loop increases with the temperature, and its logarithm reveals a characteristic 
dependence on the inverse of the square of the temperature. At higher temperatures, the 
renormalized Polyakov loop overshoots one, reaches a maximum, and then starts decreasing, 
in agreement with weak-coupling predictions. The implications of these findings are discussed. 



PACS numbers: 12.38.Gc, 11.15.Ha, ll.lO.Wx, ll.15.Pg, 12.38.Aw 



1 Introduction and motivation 

The change of state to a deconfined phase at high temperatures or densities is a very important 
phenomenon in quantum chromodynamics (QCD) and in other non-Abehan gauge theories. While 
at zero and low temperatures the physical states are color-singlet hadronic states, in the high- 
temperature limit the physical running coupling becomes small, due to asymptotic freedom, and 
one expects that the physics should be described in terms of a gas of weakly interacting quarks 
and gluons [1]: the quark-gluon plasma (QGP) [2]. These two qualitatively different phases 
should be separated by a phase transition or a crossover, which has been searched for in an 
extensive experimental heavy-ion collision programme since the 1980's. The results obtained at 
SPS, RHIC and LHC during the last decade show, indeed, convincing evidence for the creation of 
a new state of matter at temperatures about 160 MeV, which behaves as an almost ideal fluid [3]. 
The experimental research on the QCD phase diagram will be continued and extended at FAIR, 
NIC A and J-PARC. 

On the theoretical side, however, the quantitative understanding of the QCD plasma is still 
an open problem. One of the reasons for this is that the deconfined plasma retains some non- 
perturbative features even in the limit of high temperatures T. In particular, the presence of 
severe infrared divergences in weak-coupling expansions for thermal gauge theories leads to non- 
analytical properties of the perturbative series for various physical observables, and to a break- 
down of the correspondence between loop expansions and expansions in powers of the coupling [4] . 
As a consequence, the long-wavelength modes of the QGP are strongly coupled at all tempera- 
tures, and thus cannot be treated perturbatively — see ref. [5] for a review. Finally, at the typical 
temperatures probed in experiments, the physical coupling of QCD turns out to be relatively 
small, but not extremely so, and perturbative predictions fail close to the deconfinement temper- 
ature [6]. 

For these reasons, the theoretical study of the QGP at temperatures close to the deconfining 
transition is usually addressed with non-perturbative methods, including, in particular, numerical 
simulations on the lattice [7]. During the last decade, lattice computations of the equation of 
state in QCD with light dynamical quarks have reached high levels of precision, and showed that 
the deconfinement at finite temperature and vanishing quark chemical potential (for physical 
values of the quark masses) is a crossover, rather than a genuine phase transition. In fact, in 
QCD with quarks of finite mass there is no exact symmetry-breaking pattern to characterize the 
deconfinement. 

By contrast, pure SU(A^) Yang-Mills theories (which capture most of the qualitative features 
of the physics of deconfinement) provide much a cleaner theoretical setup: in the Euclidean for- 
mulation, it is easy to see that the Lagrangian of SU(A^) Yang-Mills theories at finite temperature 
is invariant under a global symmetry associated with the center of the gauge group Z^v [8]. The 
order parameter for this symmetry is the trace of the temporal Wilson line, or Polyakov loop [9]: 



In the thermodynamic limit, the ground-state expectation value of L is exactly vanishing in the 
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low-temperature phase, while it becomes non-zero above the critical deconfinement temperature 
Tc, signaling the spontaneous breakdown of center symmetry. Although L, per se, is not a 
physical observable, it can be interpreted as the trace of the propagator of an external, infinitely 
massive probe color charge located at x: a vanishing L in the ZTv-symmetric ground state at 
low temperatures T < means that the expectation value of a static color charge is zero, and 
hence the system is confined. On the contrary, L is non-zero in the high-temperature phase 
at T > Tc, corresponding to a finite free energy for the probe color charge in the deconfined 
phase. Thus, L has the meaning of an order parameter for the finite-temperature deconfinement 
transition in Yang-Mills theory. Another possible order parameter for the transition is given 
by the two-point Polyakov loop correlation function: across the phase transition, it changes 
from confining to exponentially screened. The Polyakov loop correlation function extracted from 
lattice simulations at finite temperature is often used as an input for effective potential models 
for quarkonia [10]; however, certain subtleties related to the connection between the real- and the 
imaginary-time formalism, and to the spectral decomposition into singlet and octet contributions 
to the corresponding free energies have recently been pointed out in the literature [11]. 

Note that the free energy associated with the bare Polyakov loop defined by eq. (1) is a 
divergent quantity, and hence needs to be renormalized [12]. 

In general, in SU(A^) Yang-Mills theory the Polyakov loop is an order parameter for a probe 
charge in a generic irreducible representation of the gauge group with non-zero A^-ality (i.e., 
a representation transforming non-trivially under the center of the group). The free energy 
associated with charges in different irreducible representations is expected to be proportional to 
the eigenvalue of the corresponding quadratic Casimir operator (C2). This is called "Casimir 
scaling": it is not specific to Polyakov loops, and indeed it has been studied for various other 
observables [13]. For the Polyakov loop, perturbative calculations predict Casimir scaling to hold 
at the lowest orders [14] (deviations from Casimir scaling are predicted to occur only at 0{g^)). 

In this work, we study the behavior of bare and renormalized Polyakov loops in non-Abelian 
gauge theories with a different number of colors, from 2 to 6, discussing various renormalization 
methods, and comparing our results to those of recent, similar studies for SU(3) [15-17] and 
SU(2) [18, 19] Yang-Mills theories. In particular, we investigate the features that emerge when 
N is large. The motivations for looking at the limit of a large number of colors are manifold. 
First of all, the large-A^ limit of QCD at fixed 't Hooft coupling A = g'^N and fixed number of 
flavors Nf is known to lead to dramatic mathematical simplifications [20] . For the phase diagram 
of QCD-like theories, the large- iV limit has also interesting implications for new phases at high 
density [21]. Furthermore, it plays a technically crucial role in holographic computations, inspired 
by the conjectured equivalence of maximally supersymmetric Yang-Mills theory with = 4 
supercharges in four dimensions and supersymmetric type IIB string theory in a 10-dimensional 
AdS^ X spacetime [22]. This conjecture relates the large- limit of the strongly coupled 
gauge theory to the classical gravity limit of string theory in a five-dimensional anti-de Sitter 
spacetime, which can be studied analytically. While at zero temperature the A/" = 4 theory is 
qualitatively very different from QCD, there are arguments suggesting that at finite temperature 
the two theories should share at least some qualitative (or semi-quantitative) physical features [23] . 
Calculations based on the gauge/string duality have also been extended to various other models. 



2 



which mimic the features of QCD either by breaking exphcitly some of the symmetries of the 
= 4 theory using some additional ingredients ( "top-down" approach) , or by constructing some 
ad hoc five-dimensional gravity model, which should reproduce the properties of QCD ( "bottom- 
up" approach) . These models are often used to study analytically certain features of the strongly 
coupled quark-gluon plasma [24]. 

One important technical aspect in all holographic computations is that they are based on 
the approximation of an infinite number of colors in the gauge theory: this limit allows one to 
neglect loop effects in the dual string theory, i.e. to reduce it to its classical limit. Recent lattice 
studies have showed that the large- limit is indeed a good approximation for the physical SU(3) 
case, both as it concerns spectral and thermal observables [25]; remarkably, this also holds for 
theories in 2 + 1 spacetime dimensions [26]. However, the validity of the infinite- approximation 
is, in general, a non-trivial issue, which can depend on the observable considered, and should be 
studied on a case-by-case basis. 

In the context of gauge/string duality, the behavior of the renormalized Polyakov loop as 
a function of the temperature has been recently discussed in refs. [27, 28]. In particular, in 
ref. [27] it was argued that, in strongly coupled theories with a holographic dual, the renormalized 
Polyakov loop should be monotonically increasing with T. This is in contrast with perturbative 
computations [14], which predict that the leading-order correction to the free limit is positive, and 
hence that the renormalized loop L''™ should tend to unity from above in the high-temperature 
limit. However, it should be noted that these two theoretical predictions are expected to hold 
in the strong- and in the weak-coupling regime, respectively. A holographic prediction for the 
renormalized Polyakov loop was worked out analytically in ref. [28], using a simple holographic 
model with one deformation parameter [29]. This work found that, at the leading order in a 
high-temperature expansion, the logarithm of the Polyakov loop in the strong coupling regime 
should be given by the sum of a constant plus a term proportional to (Tc/T)^, an effect which 
has also been observed and discussed in refs. [30,31]. 

The properties of renormalized Polyakov loops in theories based on different gauge groups 
are also of interest for effective models of the quark-gluon plasma in the region near Tc, see 
refs. [16,32,33] and references therein. In particular, the behavior in the large- A^ limit may reveal 
analogies with the third-order transition that one finds in 1 + 1 dimensions [34]. Moreover, at 
large A^ one expects that different irreducible representations become equivalent, up to 0{1/N) 
corrections: for example, the two-index symmetric and antisymmetric representations are ex- 
pected to be equivalent for A^ — t- oo. Furthermore, using the group theoretical tools of composite 
representations [35,36] (see the appendix A for details), it is possible to show that in the large- A^ 
limit the eigenvalue of the quadratic Casimir remains 0{N). 

Finally, the finite-temperature properties of strongly coupled gauge theories based on different 
gauge groups and with dynamical fermions in various representations are also interesting for 
extended technicolor models [37]. 

With this motivation, in this work we address a first-principle lattice study of Polyakov loops 
at finite temperature in SU(A^) gauge theories with a different number of colors A^, and for several 
irreducible representations. In particular, we consider the twelve lowest non-trivial irreducible 
representations of each gauge group, and investigate the Casimir scaling at temperatures close to 
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the deconfinement transition. Then we define non-perturbatively renormahzed Polyakov loops, 
discussing various renormahzation methods that have recently been proposed in the literature. 
While all our computations are performed in the setup of the pure Yang-Mills theory, it is worth 
remarking that, in the 't Hooft limit, the dynamics of gluons dominates, with the contributions 
from virtual quark loops suppressed by powers of the large- limit of QCD is a unitary 

quenched theory, and by virtue of this, in this limit it is legitimate to consider only the glue 
sector of the theory on the lattice. This allows one to avoid the complications arising from lattice 
fermions, and to achieve a smoother approach to the planar limit (the leading-order finite-iV 
corrections in the glue sector are proportional to 1/N'^). 

In section 2 we define the setup of our lattice computations and the method to extract the 
renormahzed Polyakov loop free energies. Our results are presented in section 3, while in section 4 
we discuss their implications, and summarize our findings. Some useful group-theoretical formulae 
are listed in the appendix A. 

Preliminary results of this study were presented in ref. [38]. 

2 Lattice simulation setup 

Our numerical simulations are based on the regularization of SU(A^) Yang-Mills theories with N = 
2, 3, 4, 5 and 6 colors on a four-dimensional Euclidean hypercubic, isotropic lattice A of spacing 
o, with periodic boundary conditions in all directions. We use natural units {h = c = = 1), so 
that the temperature equals the inverse of the size of the system in the compactified Euclidean 
time direction: T = l/{aNt), and denote the spatial volume of the lattice as V = (aNs)^. 
For most of our simulations at finite temperature, we used lattices characterized by an aspect 
ratio Ng/Nt > 4, which provides a good approximation of the thermodynamic limit [39]. The 
fundamental degrees of freedom in the lattice regularization of the theory are a discrete (and 
finite, if one considers a finite hypervolume) set of Ufj,{x) matrices in the N x N representation 
of the group, which are defined on (and represent parallel transporters along) the oriented bonds 
between nearest-neighbor sites on the lattice. The functional integral defining the continuum 
partition function of the system is traded for a well-defined, finite, multi-dimensional ordinary 
integral: 



zeA 0=1 

where dUa{x) is the Haar measure for each Ua{x) G SU(A^) link matrix, and denotes a gauge- 
invariant lattice action. The simplest choice for is given by the Wilson gauge action [40]: 




(2) 





(3) 



xeA l<fi<i/<4 '- 



with /3 = 2N/gl 



Iq and: 



U^i^) = U^{x)U,{x + a/i)f/t(x + auplix). 



(4) 
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However, in our study we used the tree-level improved gauge action [41,42]: 

'S'imp = /3 



x6A l<fj.<u<A 



--^ReTY 
2 N 



where: 



Un{x)Uu{x + afi)Uu{x + afi + ai))uUx + 2ai')Ul{x + at')Ul{i 



(5) 



(6) 



Assuming that the Uf^{x) group variables are related to the continuum gauge fields ^^(x)*^ via: 
U^{x) = exp[igQA^{x + ajl/2)ta], it is straightforward to show that both 5w and ^i^p tend to 
the Yang-Mills action in the continuum limit a — t- 0, but the tree- level improved action defined 
by eq. (5) is characterized by smaller discretization effects than those of the Wilson action. 
Expectation values of gauge-invariant physical observables O on the lattice are defined by: 
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and can be estimated numerically by Monte Carlo sampling over a finite set of {Ua{x)} configura- 
tions; in the following, we denote the number of configurations used in our computations as n^^^^f. 
The algorithm we used to generate the configurations is based on a 3 + 1 combination of local 
overrelaxation [43] and heat-bath [44] updates on N{N - l)/2 SU(2) subgroups of SU(iV) [45]. 
The parameters of our lattice simulations are shown in tab. 1. 
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Table 1: Parameters of the lattice simulations used in this work. denotes the number of colors, 
Ng and Nt are the number of sites along the space-like and time-like sizes of the lattice, is 
the number of /3-values that were simulated, in the /3niin < /3 < /?max interval. For each set of 
parameters, the number of thermalized configurations, that we used in our numerical estimates, 
is shown in the last column. 

Converting the simulation results to physical units requires a definition of the lattice scale. 
In order to set the scale for our simulations with the improved action, we calculated the T = 
static potential in lattice units from expectation values of Wilson loops {W{r,L)): 

{W{r,L-a)) 



V{r) 



~^ lim In ■ 



{W{r,L)) 



(8) 
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In particular, we extracted the potential from Wilson loops defined from smeared links, using five 
levels of smearing for the spacelike links (leaving the timelike links unsmeared). The values of 
V{r) thus obtained are then fitted to the Cornell potential: 

V{r) =ar + Vo + ^, (9) 
r 

enabling one to extract a (as well as Vq and 7) in lattice units; statistical errors are estimated 
with a jackknife analysis. All fits give Xrod values close to 1, and the 7 parameter is always very 
close to the bosonic string prediction: 7 = — 7r/12 [46] (see fig. 1 in ref. [38]). 

Note that this non-perturbative definition of the scale is not unique: in general, it would be 
equally legitimate to define the value of a (for a given /?), using the lattice results for a different 
dimensionful physical observable — for example, the critical temperature Tc [47]. On a finite- 
spacing lattice, different physical observables are generally affected by different discretization 
artifacts, and hence lead to slightly different definitions of the scale. This ambiguity is a systematic 
effect in the scale determination, but the associated relative uncertainty is numerically small, and 
vanishes in the continuum limit a — )■ 0. 

On the lattice, the trace of the bare Polyakov loop in the irreducible representation r is defined 

as: 

Nt 

Tr H Ut\x,ant) , (10) 

nt=l 

where g^^^ denotes the value of the group element g in the irreducible representation r. Note 
that the matrix elements of a generic g^"^^ can be easily obtained from those of g in the defining 
representation, by means of basic relations of representation theory. In particular, the characters 
of group elements in different irreducible representations can be easily expressed using Young 
calculus and the Weyl formula [48,49] (see the appendix A for details). 



3 Results 

3.1 Setting the scale 

To determine the scale for our simulations with the tree-level improved lattice action, we fit our 
results for a^a (as extracted from Cornell fits of the T = potential using smeared Wilson loops) 
at the largest couplings to the functional form: 

aV = exp{-[^o + ^i(/3 - /So) + A2W - Po? + MP " Pof]}, (H) 

where f3 = 2N/gQ, and /3o is an arbitrary reference value in the /3-range of our simulations. 

As an example, fitting the SU(3) data taken from ref. [42] to eq. (11) (choosing /3o = 4.3) 
yields: 

aV = exp (-2.660(12) - 3.145(66) • (/3 - 4.3) + 0.97(11) • (/3 - 4.3)^ - 0.33(26) • (/3 - 4.3)^} , 

(12) 
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with Xrod — 0.34. The corresponding data, together with the fitted curve, are shown in the top 
panel of fig. 1. 

Similarly, our data for the SU(4) gauge group yield: 

9)2 - 0.320(55)(/3 - 9)^} for /3 < 8 

for /? > 8 
(13) 

1, while for SU(5) we obtain: 

2 _ f exp (-3.021(15) - 0.682(17)(/3 - 13) +0.214(30)(/3 - 13)2} for /3 < 12.7 
" ^ ~ \ exp{-0.636(35)/3 + 5.28(45)} for /3 > 12.7 ' ^ ' 

with Xred — 3.41, see the bottom panel in fig. 1. 
3.2 Casimir scaling 

The first issue that we investigated is Casimir scaling of bare Polyakov loops, i.e., whether the free 
energy associated to bare Polyakov loops in a given irreducible representation r is proportional 
to the eigenvalue of the quadratic Casimir (C2) of that representation. To study this problem, we 
rescaled the loop free energies by the ratio of the Casimir in the given representation over the one 
in the fundamental representation /. This corresponds to raising the values of the loops to the 
power l/fi, where d = {C2)r/ {C2) f (the values of d are reported in the appendix A). Our results 
for the SU(4) gauge theory are displayed in fig. 2, which shows the values of L^/*^ for the twelve 
different representations, as obtained from simulations with the tree-level improved action on a 
lattice with Nt = 5 and Ng = 20 sites along the Euclidean time and spatial directions, respectively. 
Note that, in this plot, the data are displayed as a function of /3 = 2N/ Qq-. since the bare loops do 
not depend only on the temperature T, but also on the bare coupling qq, it is natural to display 
these values (from simulations at fixed Nt) as a function of (3. This also allows one to avoid 
introducing any potential ambiguity related to the definition of the temperature scale. If Casimir 
scaling holds, then this rescaling should make the values corresponding to higher representations 
collapse onto those of the fundamental representation (for which d = 1). 

Our results show an approximately perfect Casimir scaling in the deconfined phase, for all the 
representations that we considered. Although the bare values of loops in different representations 
vary by orders of magnitude, rescaling their free energies according to the corresponding quadratic 
Casimir eigenvalues makes them fall onto the same, universal curve. Our data show that the only 
significant deviations from this behavior (apart from the obvious ones in the confined phase, where 
Casimir scaling is not expected to hold) are visible for strongly suppressed high representations, 
which are most sensitive to finite-volume effects. For example, the rescaled bare loops in the 
representations denoted as 20", 35, 50 and 56 show significant deviations from the curve of the 
other data for temperatures T < 1.75 Tc, while they collapse on that curve at higher temperatures 
(for L^/'^ > 0.2). This is simply due to the fact that, for these representations, for T < 1.75 Tc 
the expectation value of the corresponding loops in the thermodynamic limit is smaller than the 



2 _/ exp (-3.894(38) - 1.21(14)(/3 - 9) - 0.41(16)(/3 - 
" \ exp{-l. 165(29)/? + 6.54(23)} 

with x^d ~ 1-22, and are shown in the central panel of fig. 
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Figure 1: The top panel shows a fit of the results for the string tension in lattice units in the 
SU(3) gauge theory, taken from ref. [42], to the functional form in eq. (12). The central and 
bottom panel display the fits of our results for the string tension in lattice units to eq. (13) in 
SU(4) and SU(5) Yang-Mills theories, respectively. 
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(non- vanishing) average vahie of \L\ computed on a lattice of finite vohime. This is the same 
effect that, on any finite lattice, is responsible for the non- vanishing values of \L\ in the confining 
phase. 

Fig. 3 gives evidence of this: the left panel shows our results for bare Polyakov loops in the 
fundamental representation of SU(4), obtained from lattices of two different spatial volumes, V = 
(20a)^ and V = (24a)^. The results of the two sets of simulations are compatible with each other 
in the deconfined phase (signaling that finite- volume corrections to the critical value of /3 are small 
for both ensembles), whereas the data obtained from the larger lattice are strongly suppressed in 
the confining phase, in agreement with the expectation that the average Polyakov loop is exactly 
zero in the thermodynamic limit. The right panel shows the same comparison, for loops in the 
representation of size 56: for high-dimensional representations like this, the thermodynamic limit 
value of the Polyakov loop is very small, even in large regions of the deconfined phase, and thus 
it is overwhelmed by finite-size artifacts on the lattices that we considered. As fig. 2 shows, for 
such representations it is only at very large values of (3 (i.e., at very high temperatures) that the 
contribution surviving the thermodynamic limit becomes dominant over finite-volume artifacts. 

Our results for bare Polyakov loops in different representations (rescaled by dividing the 
respective free energies by the factor d, proportional to the quadratic Casimir eigenvalue of the 
corresponding representation) for the SU(2), SU(3), SU(5) and SU(6) theories are displayed in 
fig. 4: they reveal the same behavior observed for the SU(4) gauge group. Furthermore, comparing 
the plots of the rescaled bare loops for different groups, one also observes that, when N grows, the 
numerical study of higher representations simplifies, in the sense that they tend to be less sensitive 
to finite-volume effects. This is related to the fact that, in general, for N ^ oo the quadratic 
Casimir grows only linearly with A^, and with the number of fundamental and anti-fundamental 
indices out of which a generic representation is built (see the appendix A for a discussion). For 
example, while the d factor for the highest SU(2) representation considered here (i.e., for the 
twelfth lowest, non-trivial) is equal to 56, its value for the twelfth SU(6) representation is less 
than 6. 

Note that, deep in the weak-coupling region, one could compare the simulation results with 
the predictions from lattice perturbation theory. In particular, for the Wilson action the latter 
have been known for many years [50]. However, since our simulations are based on the tree-level 
improved action, rather than the Wilson action, we did not perform such a comparison. 

3.3 Renormalized Polyakov loops 

Finally, we present our results for the renormalized Polyakov loops, restricting our attention 
to loops in the fundamental representation. Our renormalization procedure is based on the 
determination of the constant term Vq in the T = interquark potential extracted from the 
lattice, at each value of the bare coupling. More precisely, we define the renormalized Polyakov 
loop as: 



where aVo is obtained from our fits of the interquark potential. Note that, in the expression 
above, the renormalized loop L""" is expected to depend only on the physical temperature T, 
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Figure 2: Temperature dependence of bare SU(4) Polyakov loops in different representations, after 
dividing their free energies by (a quantity proportional to) the eigenvalue of the corresponding 
quadratic Casimir {C2)r- This plot displays the results we obtained from simulations with the 
tree-level improved action, on lattices with iVj = 5 and Ns = 20 sites along the compactified time 
and spatial directions, respectively. 
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Finite-volume effects on bare Polyakov loops 
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Figure 3: Left-hand side panel: Comparison of bare SU(4) Polyakov loops in the fundamental 
representation, obtained from lattices of different volumes: in the confined phase, the results 
tend to zero in the thermodynamic limit. Right-hand side panel: Loops in high-dimensional 
representations (such as the 56, displayed in this plot), whose expectation values are strongly 
suppressed, are particularly sensitive to finite-volume artifacts. 



while the bare one L depends both on and on iVj. On the other hand, Z depends only on g^. 

For SU(4), in the temperature region that we are most interested in (i.e., in the deconfined 
phase, close to Tc), our fits show that an accurate parameterization of ^(ffg) is of the form: 



Z{gl) = exp[-0.166(21)52 +0.259(28)5^] 



(16) 



for 5q < 0.8; the quoted errors are conservative. Using eq. (16) to renormalize the bare loops 
obtained from our simulations with A^^ = 5, A'^s = 20, we obtain the renormalized Polyakov loop 
values displayed in the top panel of fig. 5, in which the displayed errorbars also include an estimate 
of the systematic uncertainties related to scale setting and renormalization prescription choice 
(which are discussed below) . The inset shows a comparison of our data over a broader range (with 
extrapolation in the scale setting and in the parameterization of the renormalization constant) 
with the perturbative prediction for this gauge group, taken from ref. [14]. In particular, the 
upper solid curve is obtained using one-loop estimates for the coupling and Debye mass, whereas 
the lower dashed curve is obtained from two-loop estimates of these quantities [51]. The figure 
shows that the renormalized loop takes a value close to 1/2 for T — t- ^ and increases with the 
temperature, overshooting 1 at T ~ 2>ATc. Extrapolating our parameterizations for the scale and 
for Z to a range of small couplings values (in which we have not performed non-perturbative 
computations of the T = interquark potential) , we find that the renormalized fundamental loop 
in the SU(4) theory reaches a maximum (about 1.07) at temperatures around 30Tc, then starts 
decreasing and approaching the next-to-next-to-leading order perturbative prediction, which, for 



11 



SU(2) 



SU(3) 



— -©- ❖ 



^ ^ ^ ^ 



-V V 

♦ t * * * 



1 ! I 



6 8 10 12 14 16 18 

.1.1.1. 



1.4 1.6 1.8 2 2.2 24 2.6 2.8 3 3.2 34 3.6 



• 


fundamental 




adjoini 




4 




5 




6 




7 




8 




9 




10 




11 




12 




13 


.1.1. 


3. 


8 4 4.2 



s 


* * * * 

# • 




• 


fundamental 




■ 


sextet 






adjoint 






decuplet 






15 






15' 






21 






24 






27 






28 







35 




□ 


36 





4 4.2 4.4 4.6 4.8 5 5.2 54 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2 7.4 7.6 7.8 



SU(5) 



SU(6) 



• 


fundamental 






dectiplet 






15 






adjoint 




< 


35 






40 






45 






50 






70 






70' 




O 


75 






105 





p Kf 

12 12.5 13 13.5 14 14.5 15 15.5 16 16.5 




• 


fiindamental 




15 




20 




21 




adjoint 




56 


> 


70 




84 




105 




105' 


o 


120 




126 



17 17.5 18 18.5 19 19.5 20 20.5 21 21.5 22 22.5 23 23.5 24 24.5 25 25.5 



Figure 4: The top left panel shows results analogous to those in fig. 2, but for the SU(2) Yang- 
Mills theory. The inset shows the convergence to a universal curve, in a parameter range where 
finite-volume effects cease to dominate the results for higher representations. Similarly, the top 
right panel show the corresponding results for the SU(3) gauge group, whereas the bottom left 
and bottom right plots display the results for the SU(5) and SU(6) theories, respectively. 
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the SU(iV) Yang-Mills theory, reads [14]: 



where g denotes the physical coupling, and m^; is the Debye mass. The behavior we observe in 
our SU(4) data is consistent with the results obtained for SU(3) in previous studies [15-17]. 

Similarly, our results for the SU(5) gauge group are based on the following parameterization 
ofZ(ff2): 

Z(92) = exp[0.4115(26)5o'], (18) 

for <7o — ^-^j ^^^'^ displayed in the bottom panel of fig. 5. Similarly to the case of four colors, 
also in the SU(5) theory the renormalized loop has a value close to 1/2 for T — )• , and increases 
up to values larger than 1. 

Finally, in fig. 6, we show (minus twice) the logarithm of the renormalized fundamental 
loops for the SU(4) and SU(5) gauge groups, as a function of the inverse of the square of the 
temperature. As it was already observed in the case of the SU(3) theory [30], at temperatures 
between Tc and a few times Tc, the logarithm of the renormalized Polyakov loop appears to be 
of the form: ^ 

- 21nL'-'=" = m ( ^ ) + q. (19) 



In fig. 6, the straight lines are fits (in the temperature ranges shown in the plot legends) to eq. (19), 
which yield m = 1.1166(55), q = -0.0959(11), with xL = 0-004 for SU(4) and m = 1.4283(62), 
q = -0.3056(15), with xL = 0.003 for SU(5). The small xL values are due to the fact that 
the errorbars affecting our numerical results are dominated by the systematic uncertainties, for 
which we could only provide a crude, but conservative, estimate. 

Note, however, that the statement, that the logarithm of the renormalized Polyakov loop is of 
the form appearing on the right-hand side of eq. (19), is a scheme-dependent one. For example, 
redefining the renormalized Polyakov loop free energy with the addition of a constant, would 
introduce an additive contribution 0{T~^) to the logarithm of the renormalized loop. We find 
that the statement holds for the renormalization scheme that we discussed here (see also ref. [15]). 

In view of this observation, one may wonder, whether there are arguments supporting our 
scheme choice, rather than others. As discussed above, our renormalization scheme for the 
Polyakov loop is based on the subtraction of the constant term appearing in the T = po- 
tential between two static sources. This reduces the form of the renormalized confining potential 
to: 

V(r) =ar + - + 0(r~^). (20) 
r 

The functional form appearing on the rigt-hand side of eq. (20) can be derived (at the leading 
order in an expansion around the large-distance limit) from an effective bosonic string model 
for confinement [46]. Various recent works (see, e.g., ref. [52] and references therein) show that 
Lorentz-Poincare symmetries constrain the first few terms in the expansion of the effective string 
action to equal those that are obtained expanding the Nambu-Goto action [53] , while corrections 
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Figure 6: Similarly to what was observed for the SU(3) theory [30], also the logarithm of the 
SU(4) (left-hand side panel) and SU(5) (right-hand side panel) renormalized Polyakov loops 
exhibits a characteristic r~^-dependence in the deconfined phase, up to temperatures of a few 
times Tc- Note that the errorbars include conservative estimates of the systematic uncertainties 
(see subsection 3.4). 



only appear at high orders in 1/r. The fact that the Nambu-Goto string provides a good effective 
model for the confining potential is also confirmed by extensive numerical evidence from lattice 
simulations, both for SU(iV) Yang-Mills theories [54] and for theories based on smaller gauge 
groups [55]. Following ref. [56], it is then natural to define a renormalization scheme yielding 
a T = interquark potential with a vanishing constant term, eq. (20), and to apply it to the 
renormalization of the Polyakov loop. 

In ref. [28], a holographic prediction for the renormalized Polyakov loop was computed, using 
a model with one deformation parameter [29]. The result reads: 



L^™(T) = 6iexp 



Vtt— Erfi I — I - exp 



(21) 



where h\ and 62 are two coefficients that can be fitted, and Erfi denotes the imaginary error func- 
tion. At the leading order in a high-temperature expansion, eq. (21) predicts that the logarithm 
of the Polyakov loop would be given by the sum of a constant plus a (Tc/T)^ term, as observed 
in the numerical data. 

In the literature, it was suggested that the dependence of In L'"''^ on could be due to a 
non-perturbative contribution from a gluon condensate [30,31]. Similar arguments have been 
invoked to explain the behavior of the interaction measure A at temperatures of the order of a 
few times Tc [57]: in all SU(iV) gauge theories, both in D = 3 + 1 [25] and in L> = 2 + 1 [26] 
spacetime dimensions, A appears to be proportional to T^. 
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3.4 Systematic uncertainties 

Apart from the precision limits related to the finiteness of our statistical samples, the main 
systematic uncertainties affecting our study include: ambiguities in the scale determination, 
renormalization prescription dependence, effects due to the volume finiteness, and finite-cutoff 
effects. Let us discuss each of them in turn. 

In the temperature range of interest in this study, a reliable definition of the temperature scale 
is necessarily non-perturbative, and — as discussed above — requires the choice of a dimensionful 
physical observable of reference. As different observables are generally affected by different lattice 
artifacts, this leads to slight ambiguities in the definition of the scale; however, this systematic 
effect becomes negligible at small lattice spacings. A potentially more severe ambiguity is re- 
lated to the functional form that one can choose to parameterize the data to be fitted. Rather 
than interpolating our simulation results with arbitrary, arbitrarily complicated functions, in the 
present study we tried to use physically motivated functional forms, with a minimal number of 
parameters, and estimated the systematic uncertainty related to scale setting by comparing the 
results with different parametrizations, at various values of the lattice gauge coupling. 

A potentially large systematic ambiguity in our computation is related to the choice of the 
Polyakov loop renormalization method. In the present work, we followed the approach already 
used in a similar study for the SU(2) gauge theory [18]. Other related studies discuss different 
renormalization methods, which lead to roughly compatible results. In particular, the authors of 
ref. [17] discussed a comparison of a renormalization method based on the QQ potential (similar 
to our prescription) with an iterative renormalization (based on simulations on lattices of different 
spacing and at the same temperatures): these two methods appear to be compatible with each 
other, although the latter has the drawback of leading to an accumulation of statistical errors, 
particularly at temperatures close to Tc- A different renormalization method was suggested in 
ref. [16]: there, the idea is to extract the free energy of the renormalized Polyakov loop at a 
given temperature T, by identifying the iVf-independent contributions to the free energy F of 
bare loops extracted from simulations on lattices of different spacing: 

F = NtF"'^ + + F'^^'/Nt + ..., (22) 

where F"^'" and F'''* respectively denote the coefficient of the contribution to the bare free energy 
that diverges in the continuum limit, and the coefficient of the leading term due to lattice artifacts. 
A problem with this method, however, is that, in order to keep the temperature T = l/{aNt) fixed, 
the lattice spacing a is obviously different for each simulation at a different Nf. Since a is tuned 
by varying the bare coupling go, this implies that F'*'", F"""" and F'^'*, which generically depend 
on go, are not held fixed when T is fixed. Yet another renormalization method was proposed in 
ref. [19], following the fixed-scale approach [58]: the idea is to fix Z at only one value of the bare 
coupling go, and then to vary the temperature in the lattice simulations by varying iVj at fixed 
spacing a (i.e., at fixed bare coupling). A potential drawback of this method, however, is that it 
does not allow one to vary the temperature continuously. Aspects related to the renormalization 
of Wilson lines have also been discussed in ref. [56]. To get a rough estimate of the systematic 
uncertainty associated with the choice of a renormalization method, we compared the difference 
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between various methods, at different temperatures, both in our data and in the results available 
in the literature. 

Finally, finite-volume and finite-cutoff effects appear to be under control in our study. In 
particular, the results of our simulations show that, for > 4, deviations from the thermody- 
namic limit in the deconfined phase are clearly visible only for high representations, whereas they 
appear to be negligible for the fundamental representation (see fig. 3). In fact, the lattices used 
in the present study are characterized by an aspect ratio Ng/Nt > 4, which is known to provide 
a good approximation of the thermodynamic limit in the temperature range of interest [39]. As 
for finite-cutoff effects, unfortunately we could not repeat all of our calculations on finer lattices, 
hence we are unable to perform a continuum extrapolation of our results. The systematic error 
due to cutoff effects on lattices with Nt = 5, however, is expected to be rather small for simu- 
lations with the improved action that we used: previous studies for the SU(3) gauge group [17] 
showed no significant discrepancies between Nt = 4 and Nt = 8. 

Adding up the various sources of systematic uncertainties in quadrature, the total relative 
errors on our renormalized Polyakov loops are in the range between 1% and 5% for SU(4), and 
between 1.5% and 8% for SU(5). 

4 Discussion and conclusions 

Our main findings can be summarized as follows. 

1. For all gauge groups, and for all the representations considered in this work, the bare 
Polyakov loops show excellent Casimir scaling, for all values of the coupling (or, equivalently, 
of the temperature), down to the deconfining transition. The only deviations, that our data 
reveal, can be explained in terms of finite-volume artifacts: they especially affect the high 
representations, particularly close to Tc (while they become negligible at sufficiently high 
temperatures), and can be reduced by increasing the spatial volume of the system. As we 
mentioned, in the literature, several works have reported evidence for Casimir scaling in 
SU(A^) Yang-Mills theories, including, in particular, for the T = string tension associated 
to the potential between two static sources in a given representation [13]. This observation 
has an interesting implication related to the large- A'^ limit. As it is well-known, for pure 
SU(A'^) Yang-Mills theories, the expansion around the N ^ oo limit can be organized in 
a series of powers of 1/A^^, i.e. it does not contain odd powers of 1/A^. As discussed 
in ref. [59], this expectation seems to be at odds with the numerical evidence of Casimir 
scaling of fc-string tension from lattice simulations [13], since, in general, if Casimir scaling 
holds, then the leading finite- corrections are 0{1/N). The resolution of this apparent 
paradox, however, was recently pointed out in ref. [60], and is based on a cancellation of 
terms involving odd 1/N powers in the spectrum of open string states. Similar arguments 
were also discussed in ref. [61]. 

2. In the thermodynamic limit the renormalized fundamental Polyakov loop is vanishing in 
the confined phase, and jumps to a finite value at the critical temperature, compatibly 
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with the first-order nature of the deconfining transition. The hmit of L''°" for T — )• is 
a number close to 1/2, which, interestingly, is the value that is obtained analytically for 
— >• oo in 1 + 1 spacetime dimensions [34]. For T > T^, the renormalized Polyakov loop is 
at first growing with the temperature (in the regime in which the plasma is most strongly 
coupled), it overshoots the value 1 at temperatures around ST^, reaches a maximum, and 
then eventually starts decreasing, in agreement with the perturbative predictions [6]. 

3. In the deconfined phase, for temperatures up to approximately 3Tc or 4Tc, the logarithm of 
the renormalized Polaykov loop (in the renormalization scheme that we considered here) is 
described well by the sum of a term inversely proportional to the square of the temperature, 
plus a constant. 

4. The finite-temperature behavior of gauge theories based on different SU(A^) gauge groups 
appears to be qualitatively and quantitatively very similar (confirming previous studies 
both in 3 + 1 [25] and in 2 + 1 dimensions [26]). The precision and accuracy limits in this 
study do not allow us to extract a reliable estimate of the (small) differences between the 
various groups. 

In conclusion, our study shows that, in the deconfined phase, the Polyakov loop satisfies 
Casimir scaling, and is only mildly dependent on the number of colors A^. The independence on 
the rank of the gauge group (which has also been observed for the equation of state per gluon 
d.o.f. [25]) supports analytical approaches based on the large- limit, including, in particular, 
holographic computations. Our results for the renormalized Polyakov loop show that this quantity 
interpolates between a regime (possibly dominated by contributions of non-perturbative nature) 
in which it is increasing with T, and one in which it tends to the perturbative prediction, and 
decreases with the temperature, approaching 1 from above in the weak-coupling limit for T — )■ oo. 

For the future, we plan to extend the present study of large- gauge theories at finite temper- 
ature, by looking at other observables, which could potentially reveal a stronger dependence on 
the rank of the gauge group. Of particular phenomenological interest are transport and diffusion 
coefficients — see, e.g., ref. [62] for a review. 
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A Irreducible representations of the algebra of generators of 

SU(7V) 

In the following, we discuss the classification of irreducible representations of the algebra of 
generators of a generic special unitary group of degree N. Further details can be found, e.g., in 
ref. [63]. 

A generic irreducible representation of the algebra of generators of SU(A^) can be labelled by 
— 1 non- negative integers Ai, A2, A3, . . . , Aat-i, with: 

Ai > A2 > A3 > ••• > Aat^i > 0. (A.l) 

The [Ai, A2, A3, . . . , Aat-i] sequence can be uniquely associated to a Young diagram with rows of 
lengths Ai, . . . , A^v-i. An alternative way to identify an irreducible representation is in terms 
of its canonical label (mi, m2, . . . , m7v_i), where the mj's represent the differences in lengths 
of subsequent rows in the corresponding Young diagram: rui = Aj+i — Aj for i < N — 1, and 

Particularly interesting irreducible representations of SU(A^) include the fundamental one 
[1,0, 0, . . . , 0], of dimension A^, the trivial one [0, 0, 0, . . . , 0] of dimension 1, and the adjoint one 
[2, 1, 1, . . . , 1], of dimension A^^ _ i_ 

More in general, the dimension of an irreducible representation is given by the formula: 

^-1 ^ /._/. 

n n irr^' k = X^ + N-^, l^ = N-i, (A.2) 
1=1 j=i+i * i 

with Atv = or, equivalently: 

N-l N-li+l-l N-1 

AT n n E ("^'^ + ^ith^ j^N=ii (ti). (A.3) 

^ 1=1 i=i k=i t=i 

A common way to denote irreducible representations is via their dimension; note, however, that 
this may be ambiguous (except for SU(2)), since in general there can be inequivalent irreducible 
representations of the same dimension. For example, SU(4) has three inequivalent irreducible 
representations of dimension 20, which can be denoted as 20, 20' and 20". In such cases, our 
convention is to use the notation with the least primes for the representation with the smallest Aj 
for the minimum value of i (namely, for the representation described by a Young diagram with 
the smallest number of boxes in the top row, or in the highest row which is different from the 
other representations of the same dimension). So, for instance, for SU(3) the [3, 1] irreducible 
representation is denoted as 15, while the [4,0] will be denoted as 15'; for SU(4), the [2, 1,0] is 
denoted as 20, the [2, 2, 0] is denoted as 20', and the [3, 0] is denoted as 20". 

The A^-ality of an SU(A^) representation defines its transformation properties under the center 
of the group, Z^v, and it is given by the total number of boxes appearing in the Young diagram, 
modulo A^. Representations of vanishing A^-ality (such as the trivial representation and the 
adjoint one) are blind to the action of the transformations in the group center. 
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Given an irreducible representation r = [Ai, A2, • • • , Atv-i], its conjugate representation is: 
f = [Ai, Ai — Xn-i, Ai — AAr_2, . . . , Ai — A2], so that its Young diagram is obtained by fitting the 
diagram of the representation r in a rectangle of N rows and Ai columns, removing all the boxes 
belonging to the Young diagram of r, and turning the diagram with the remaining boxes by an 
angle vr. 

Obviously, two mutually conjugated representations have the same dimension, and, given 
that their respective characters are obtained from each other by complex conjugation, we only 
include one of them in our lists of irreducible representations. It is most natural to use the 
"barred" notation for the representation with the Young diagram with more boxes, so that, for 
example the [1, 0] representation of SU(3) is denoted as 3, while its conjugate representation [1, 1] 
is denoted as 3. 

Representations which are self-conjugate have real characters; in particular, this is always 
the case for the trivial and for the adjoint representations. Also, note that, for a self-conjugate 
irreducible representation, the canonical label is a palindrome. 

In order to discuss the large-A^ scaling of the size and quadratic Casimir of an irreducible 
representation r, it is convenient to introduce the non-negative integers / and m, which represent 
the minimum number of fundamental and anti-fundamental factors from which the representation 
r can be constructed (by tensor products). I and m can be easily obtained from the Young diagram 
of r: I IS given by the sum of the number of boxes in all columns of length not larger than N/2, 
while m is given by the sum of the number of missing boxes in all columns of length larger than 
N/2. The A^-ality of a representation is given by (/ — m) modulo N. In the large- limit, it 
is possible to show [35, 36] that characters of different representations only depend on I and m, 
and that, although the dimension of the representation r grows like A^'+™, the eigenvalue of the 
quadratic Casimir is linear in N: 

{C2)r = ^[l + m + 0(1/A)] . (A.4) 

For SU(2), all irreducible representations are self-conjugate. The Young diagram of a generic 
irreducible representation of spin j = n/2 consists of one horizontal row of n boxes; bosonic 
representations correspond to even values of n, and have vanishing A-ality, while fermionic rep- 
resentations correspond to odd values of n, and their A-ality is 1. The associated canonical label 
is (n) (with / = n, m = 0), the dimension is n + 1, the eigenvalue of the quadratic Casimir (defined 
according to our conventions) is (C2) = ?i(n + 2)/4, and its ratio with respect to the fundamental 
representation is d = n{n + 2)/3. 

For larger SU(A) groups (upt to A^ = 8), the lowest irreducible representations are listed in 
tables 2-7. 

Generically, the eigenvalues of a group element g in the fundamental representation of SU(A) 
lie on the unit circle in the complex plane, and their product is 1: 

N 

5/ = [/•diag(e^"i,e^"2,e'"3,...,e'"'^) -[/^ with: ^ = mod 27r. (A.5) 

1=1 
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Young diagram 



□ 

p 



rrm 
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iV-ality 



canonical label 



(1,0 
(2,0 

(1,1 
(3,0 

(2,1 
(4,0 
(5,0 

(3,1 
(2,2 
(6,0 

(4,1 
(7,0 



dimension 



3 
6 
8 

10 
15 
15' 
21 

24 

27 
28 
35 
36 



notes 



fundamental 
adjoint 



self-conjugate 



(C2 



4/3 
10/3 

3 

6 

16/3 
28/3 
40/3 

25/3 

8 
18 
12 
70/3 



d 



1 

5/2 
9/4 
9/2 

4 

7 

10 

25/4 
6 

27/2 
9 

35/2 



Table 2: The irreducible representations of the SU(3) gauge group studied in this work. For this 
group, the integers I and m of each representation are respectively equal to the first and second 
index in the canonical label. 



Knowing the eigenvalues of gj, it is possible to calculate explicitly the character of g in any 
irreducible representation r = [Ai, A2, . . . , Aat.i] by means of the Weyl formula [48,49]: 

detF(A) , 

Trgr = (A.6) 

detF(O) ^ ' 

where -F(A) is an x matrix with entries defined as: Fki{X) = exp [i {N + A/ — q^], with 
Aat = 0, and e*"^, e*"^, . . . e*"^ are the eigenvalues of g in the fundamental representation. 

In many cases, however, the characters in high-dimensional irreducible representations can 
be more expediently calculated, using the laws of representation composition encoded in Young 
calculus, and using the well-known fact that the character in a representation which can be 
expressed as the tensor sum (product) of two representations is equal to the sum (product) of 
the characters in the summand (factor) representations. 



A.l Casimir operators 

A Casimir operator of a Lie algebra 5 is a homogeneous polynomial of order p, lying in the 
enveloping algebra of g, T{g), and commuting with all elements of g. Given a Casimir operator 
Cp, any product of it by an arbitrary scalar factor aCp, as well as any integer power of it Cp, 
are also Casimir operators; however, the number of independent Casimir operators of a given 
algebra g is equal to the rank / of the algebra. In particular, the algebra of generators of the 
special unitary group SU(A'') has A^ — 1 independent Casimir operators C2, C3, ...Cjy, whose 
eigenvalues (Cp) can be used to classify the irreducible representations of the algebra. 
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Table 3: Same as in table 2 (with the addition of the I and m indices), but for the SU(4) gauge 
group. 



Explicit expressions for the CpS can be obtained as follows. Starting from a basis {Eij}^ .^^ ^ 
of generators of U(A^): 

[Ea,hi -E-c.d] = h,cEa,d — Sa,dEc^b, (A. 7) 

introduce a basis for the algebra of generators, denoted as I-Bij l (where both i and j run from 
1 to A^, but the element En,n element is not defined), through: 

I Ei,i - jf 2_,k=l Ek,k II i - J. 
For the generators of SU(A^), the Casimir operator of order p can then be defined as: 

1 ^ , , 

Note that, by construction, the linear Casimir operator Ci is identically vanishing on the algebra 
of generators of SU(A^), as they are all traceless. 

The eigenvalue of Cp in the generic irreducible representation labelled by [Ai, A2, . . . Aat-i] can 
be obtained in the following way (taking Aat = 0) [64]: 

1. define A = Yl^=i ^i'^ 

2. define = Aj - \/N for alH = 1, 2, . . . N; 



22 



Young diagram 



iV-ality 



canonical label 



dimension 



m 



notes 



(C2 



d 



□ 

B 



rrm 



1,0,0,0 
0,1,0,0 
2,0,0,0 

1,0,0,1 

3,0,0,0 
1,1,0,0 

1,0,1,0 

0,2,0,0 

2,0,0,1 

4,0,0,0 
0,1,1,0 

2,1,0,0 



5 

10 
15 

24 

35 
40 

45 

50 

70 

70' 
75 

105 



fundamental 



adjoint 



self-conjugate 



12/5 
18/5 
28/5 



48/5 
33/5 

32/5 

42/5 

42/5 

72/5 
8 

52/5 



1 

3/2 
7/3 

25/12 

4 
11/4 

8/3 

7/2 

7/2 
6 

10/3 
13/3 



Table 4: Same as in table 3, but for the SU(5) gauge group. 

3. define pi = N — i and li = rrii + pi for all i = 1, 2, ... A^; 

4. for all k > 2, construct the quantities: = YliLi {^i ~ Pi)'i 

5. for all A: > 2, define the coefficients: = Ylj=i 'jt^^zj^^Sj; 

6. construct the function: ip{z) = ^^^2'^^-^^' 

7. calculate the Bp coefficients from the following Taylor expansion around z = 0: 

1 - exp [-ip{z)] 



p=0 



(note that Bq = 0); 
i. compute the eigenvalue of Cp from the formula: 



(Cp) 



Bp - NB, 
V 



(A.IO) 



(A.ll) 
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Youn 


y 


diagram 


iV-ality 


canonical label 


dimension 


/ 


m 


notes 


{C2) 


d 




□ 


1 


(1,0,0,0,0) 


6 


1 





fundamental 


35/12 


1 








2 


(0, 1,0,0,0) 


15 


2 







14/3 


8/5 








3 


(0,0,1,0,0) 


20 


3 





self-conjugate 


21/4 


9/5 






□ 


2 


(2,0,0,0,0) 


21 


2 







20/3 


16/7 




: 







(1,0,0,0,1) 


35 


1 


1 


adjoint 


6 


72/35 


L 


1 


1 1 
J 1 


3 


(3,0,0,0,0) 


56 


3 







45/4 


27/7 






□ 


3 


(1, 1,0,0,0) 


70 


3 







33/4 


99/35 






□ 


5 


(1,0,0,1,0) 


84 


1 


2 




95/12 


19/7 






□ 


4 


(1,0,1,0,0) 


105 


4 







26/3 


104/35 










4 


(0,2,0,0,0) 


105' 


4 







32/3 


128/35 








1 1 


1 


(2,0,0,0,1) 


120 


2 


1 




119/12 


17/5 


L 






_u 


4 


(4,0,0,0,0) 


126 


4 







50/3 


40/7 



Table 5: Same as in table 3, but for the SU(6) gauge group. 
This gives, in particular, the following relations: 



1 
4 



^3 + ( 2 - ^ ) ^2 



54 + (2 - iV) 53 + ( 2 - -iV ) S2 



(Ce) 



10 



S5+[^-N]S^+{=^-2N]S3+{--2N]S2--Si 



56 + (3 - TV) ^5 + ( 5 - -iV ) 54 + ( 5 



N 



(A.12) 
(A.13) 
(A.14) 
(A.15) 



}^n)ss+(3-^N]S2 



-S2Ss+[--l]sl 



(A.16) 
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Young diagram 



iV-ality 



canonical label 


dimension 


/ 


m 


(1,0,0,0,0,0) 


7 


1 





(0,1,0,0,0,0) 


21 


2 





(2,0,0,0,0,0) 


28 


2 





(0,0,1,0,0,0) 


35 


3 





(1,0,0,0,0,1) 


48 


1 


1 


(3,0,0,0,0,0) 


84 


3 





(1,1,0,0,0,0) 


112 


3 





(1,0,0,0,1,0) 


140 


1 


2 


(2,0,0,0,0,1) 


189 


2 


1 


(0,2,0,0,0,0) 


196 


4 





(1,0,1,0,0,0) 


210 


4 





(4,0,0,0,0,0) 


210' 


4 






notes 



{C2) 



d 



□ 

B 

m 



cm 



rrm 



fundamental 



adjoint 



24/7 
40/7 
54/7 

48/7 



7 

90/7 
69/7 

66/7 

80/7 

90/7 
76/7 
132/7 



1 

5/3 
9/4 



49/24 

15/4 
23/8 

11/4 

10/3 

15/4 
19/6 
11/2 



Table 6: Same as in table 3, but for the SU(7) gauge group. 



S,+ [1-n\Sq + {7- 3iV) ^5 + j"^ - 5iV ) 54 + (7 - 5iV) ^3 



+ [--^N\S2-SiS2-'-Sl+(-^ + N]S:,S2+(-'^ + \N] Sl 



,(A.17) 



^8 + (4 - N) ^7 + (f - 56 + (14 - 7N) S, + (l4 - ^iV ) ^4 

+ I f - 7iv) S3 + f 4 - ^-n] S2 - S5S2 - 5453 + (-4 + N) S^S2 



+ 



11 . 25 
'2 



+ ] Si + ^Sl 



(A.18) 



In turn, the equations above lead to the following expressions for the quadratic Casimir 
eigenvalues (C2): 



(C2) = -Ai (Ai + 2) forSU(2), 



(A.19) 
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Youn 


y 


diagram 


iV-ality 


canonical label 


dimension 






notes 


{C2) 


d 




□ 


1 


fl 0) 


8 


1 





fundamental 


63/16 


1 




a 


2 


(0,1,0,0,0,0,0) 


28 


2 







27/4 


12/7 


m 


2 


(2,0,0,0,0,0,0) 




2 







35/4 


20/y 








3 


(0,0,1,0,0,0,0) 


56 


3 







135/16 


15/7 






□ 





(1,0,0,0,0,0,1) 


63 


1 


1 


adjoint 


8 


128/63 








4 


(0,0,0,1,0,0,0) 


70 


4 





self-conjugate 


9 


16/7 


1 




1 


3 


(3,0,0,0,0,0,0) 


120 


3 







231/16 


11/3 


[ 


3 


□ 


3 


(1,1,0,0,0,0,0) 


168 


3 







183/16 


61/21 






□ 


1 


(1,0,0,0,0, 1,0) 


Jib 


i 


2 




i / 5/ lb 


OK /Ci 

2b/ V 








1 1 


1 


(2,0,0,0,0,0,1) 


280 


2 


1 




207/16 


23/7 


1 


1 


1 


1 1 


4 


(4,0,0,0,0,0,0) 


330 


4 







21 


16/3 








4 


(0,2,0,0,0,0,0) 


336 


4 







15 


80/21 



Table 7: Same as in table 3, but for the SU(8) gauge group. 



(C^2) = - (A? + 3Ai - A1A2 + A2) forSU(3), (A.20) 

(C2) = I [3X1 + A2 (4 + 3A2) - 2A3 (A2 + 2) + 3Ai - 2Ai (A2 + A3 - 6)1 for SU(4), (A.21) 

o 

(C2) = ^ [2 (A? + + + Xl) - (5 + A3) A4 + A2 (5 - A3 - A4) 

+Ai (10 - A2 - A3 - A4)] for SU(5) , (A.22) 

(^2) = 3^ [5 (A? + Xl + Xl + Xl + Xl)+6 (A3 - A4) - 2A3A4 - 2 (9 + A3 + A4) A5 

+2A2 (9 - A3 - A4 - A5) + 2Ai (15 - A2 - A3 - A4 - A5)] forSU(6), (A.23) 



26 



{C2) 



1 



[3 (A? + Xl + Xl + Xl + Xl + Xl)+7 (3Ai + 2X2 + A3 - A5 - 2A6) - A4A5 

A3 (A4 + As) - (A3 + A4 + A5) Ae - (A3 + A4 + As + As) A2 
Ai (A2 + A3 + A4 + As + Ae)] forSU(7), ( 



7 



(A.24) 



{C2) = [7 (A? + Xl + Xl + Xl + Xl + Xl + A?) + 24A3 + 8A4 - 2A3A4 - 8As - 2A3AS 
-2A4AS - 24A6 - 2A3A6 - 2A4A6 - 2AsA6 - 2 (20 + A3 + A4 + As + Ag) A7 
-2A2 (-20 + A3 + A4 + As + Ae + A7) 



Note that the Casimir operators are defined up to a multiphcative constant; with the conven- 
tions fixed by the construction above, the eigenvalue of the SU(A^) quadratic Casimir operator 
in the fundamental representation is (A^^ — l)/(2iV), while in the adjoint representation it is A^. 
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